load("pairData.Rda")


##############################################
## unidimensional pairwise comparison model
##############################################
library(MCMCpack)
set.seed(580223)
out.1d.alpha <- MCMCpaircompare(pairData, 
                                theta.constraints=list('1015' = 0.25,
                                                       '1014' = "-"),
                                alpha.fixed=FALSE,
                                a=0, A=1,
                                verbose=1000,
                                burnin=20000, mcmc=100000, thin=10,
                                store.theta=TRUE, store.alpha=TRUE,
                                seed = 58203 )

save(out.1d.alpha, file="out-1d-alpha-nofraud-balanced-newIDs.Rda")

